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Abstract. We present new brightness monitoring observations of the 4 components of gravitationally lensed system 
Q2237+0305, which show detection of an intrinsic quasar brightness fluctuation at a time of subdued microlensing activity, 
between June 27 and October 12, 2003. These data were used to determine the time delays between the arrivals of the four 
images. The measured delays are t ba ss -6, tqa ~ 35, and t da ss 2 hours, so they confirm that the long history of brightness 
monitoring has produced significant detection of microlensing. However the error bars associated with the delays, of order 2 
days, are too large to discriminate between competing macro-imaging models. Moreover, our simulations show that for the 
amplitude of this intrinsic fluctuation and for photometric errors intrinsic to optical monitoring from our 1.5-m telescope or 
from the OGLE monitoring, a daily sampled brightness record cannot produce reliable lags for model discrimination. We use 
our simulations to devise a strategy for future delay determination with optical data. Nevertheless, we regard these first esti- 
mates to be significant, since they are the first direct measurements of time delays made for this system from ground-based 
observations in the visual wavelengths. Our results provide the most convincing confirmation of the gravitational-lens nature of 
Q2237+0305, and give observational justification to the extensive literature which attributes the quasar's previously observed 
brightness fluctuations to microlensing. 
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£^ | 1. Introduction of the Q2237+0305 components did not reveal any region of 

' , • • „ , , correlated brightness variations, which could be attributed to 

S/* The quadruple Q2237+0305 gravitationally lensed quasar sys- t , " . , . , „_ „ ... 

rS f ■< o j ~i j tne q Uasar intrinsic brightness fluctuation, with one exception; 

$H ' tem was among the first discovered (Huchra et al. 1985 Yee . , i.»»,i . , , , . . c 

f3 ° . , 1 1 0stensen et al.(1996) noted a synchronous brightening of all 

^ 1988 1, and has frequently been cited as an important system in „ ^ . , „„. . ^ , „ . . . . 

, ' . , . , •,, , , ,,, , four components in 1994, - in the very end of their monitoring 

which microlensing by stars in the lens galaxy should be read- ... ,, XT ,. „ . . „ . r . . . , 

° J a j program with the Nordic Optical Telescope. It was interpreted 

lly observed with events having duration of a few hundred days. ... . , . . . . ... ,. „ nnn£ 

J ° J as a possible intrinsic variation of the source, reaching 0.0006 

However microlensing studies have ordinarily required deter- .. , . ™ c , , . a- . . c 

b magnitudes per day. The authors found this to be insufficient for 

mination of time delays to permit discrimination between in- , . . .. , , , . . , , . .. 

J r determining the time delays, but noted a rather high correlation 

trinsic quasar fluctuations and mi crolensing-induced brightness , . . . . . . , , 

M & & between the components. A single measurement of time delay 

fluctuations. In the case of 02237+0305, models of the macro- „ ^^->n rvorv^ £ j i. r. • ^ i £ 

_ , .rrrrmnT.- f° r Q2237+0305 exists so far, made by Dai et al. ( 2003 ) from 

imaging from the earliest times (e.g. Schneider et al. 1988 Rix 



et al. 119921 Wambsganss & Paczynski I1994L as well as the 
most recent and most popular model by Schmidt, Webster & 



the Chandra X-ray observations. They reported a 2.7 hours time 

delay between images A and B, with image A leading. 

. When the first Q2237 microlensing fluctuations were an- 

Lewis ( 1998 1 all suggested that the time delays should be short, , . lr . n u T ^ ^ ^ nnncn j 

, — ' , , , ...... nounced in 1989 by Irwin et al. (1989), the Q0957 time de- 

of order a few hours up to 1 day, but this prediction has never . . , , , . , , „ .... „ .„ mrsn 

, „ ,, , . ' . . f . , , . , lay had already been measured by Schild & Chorfin (19861 

been confirmed by observations. This is because of the high , „ , . , . , „ 

. ... ... r ... „ and conhrmed by Vanderriest et al. (19891. Both authors re- 

optical density to microlensing for all the four quasar com- . , ., „ . , . . , . 

r , • , • ported evidence of microlensing in the Q0957 system (Gneger, 

ponents that their brightness changes are unceasing and sig- „ , _ c , , i lr , , i , . , . i lr , or ,. , . lU 

, , , _„ , , ., , , , ° Kayser, and Refsdal 1986; Vanderriest et al. 1989), but the 
mficantly uncorrected. Till recently, the available light curves . . . ^ . , . , , 
J ° claim was not completely accepted because the time delay was 



Send offprint requests to: V.Vakulik, e-mail: challenged by Lehar et al. ( 1 992) from radio data and by Press 
vakulik@astron.kharkov.ua et al. ( 1992) from re-analysis of the optical data. However re- 



2 



V.Vakulik et al.: Time delays in Q2237 



determination of the Q0957 time delay by Schild & Thomson 
fi997l . Pelt et al. fT996l and Kundic et al. fT997l firmly estab- 
lished the original Schild & Cholfin ( 1986 1 time delay and, thus 
also, the first discovered microlensing event, evidently caused 
by a star in the lens galaxy (Schild & Smith 1991). 

During the same time period, numerous events were de- 
tected in the Q2237 system and discussed in the context of mi- 
crolensing even though the time delay had not been measured, 
but was presumed to be short, (Corrigan ll991l Wambsganss, 
Paczynski, and Schneider 1990 1. Thus it would be important to 
observationally confirm this prediction in order to validate the 
microlensing inferences. However because the microlensing of 
the individual images is seen nearly continuously, it has been 
difficult to date to detect an intrinsic quasar fluctuation from 
which the time delays could be determined. A similar conun- 
drum has been discussed in the case of Q0957 by Colley and 
Schild ( 2003 1; because of the microlensing it has been difficult 
to determine an accurate time delay, but without an accurate 
time delay it has been difficult to study the microlensing. 

We report below the observation of an unambiguous Q2237 
intrinsic quasar fluctuation seen in all four images. In section 2 
we report our observations with the 1 .5-m Maidanak telescope, 
and their reduction to form the four time series records. In sec- 
tion 3 we determine the delays between the arrivals of the four 
images by two methods, and in the final sections we discuss our 
results in the context of intrinsic quasar fluctuations observed 
in other quasars, and in the context of the microlensing and 
macro-lensing models. 

2. Observations and data reduction 

CCD images of the Q2237+0305 gravitationally lensed quasar 
were taken with the Maidanak 1.5-m telescope in the frame- 
work of an international cooperative program of monitoring 
gravitational lenses, launched in 1997, (Dudinov et al. 2000b 
Vakulik et al. 2004 1. An LN-cooled BroCam CCD camera with 
a SITe ST 005A chip was available in the f/8 focal plane, giv- 
ing a scale of 0.26"/pix. Magnitude measurements of the four 
quasar components with the corresponding uncertainties are 
shown in Table where the Julian dates and the seeing con- 
ditions (FWHM values) are also indicated. 

Because of some problems with the telescope tracking sys- 
tem, we used rather short exposures, which do not exceed 3 
min. Therefore, to provide the highest photometric accuracy, 
the Q2237 images were usually taken in series, consisting of 
4 to 12 frames. The frames were averaged before being sub- 
jected to photometric processing, thus ensuring the accuracy of 
about 0.008" 1 for A, 0.018'" for B, 0.016'" for C, and 0.020'" for 
the D component. These are the mean errors, averaged over the 
entire dataset, while the error values presented in Table ^for 
every date were obtained from a comparison of the photometry 
of individual frames in the series, and may be regarded as an 
adequate estimate of the random error inherent in a particular 
quasar image component series. 

The algorithm of photometric image processing is de- 
scribed in great details in (Vakulik et al., I2004> . where the 
absence of seeing-dependent systematic photometric errors for 
images with PSF values up to 1 .4" is claimed. It is interesting 
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Fig. 1. Synchronous variations of brightnesses of the A-D im- 
ages observed in Q2237+0305 between June 27 and October 
12, 2003 with the Maidanak 1.5-m telescope. R magnitudes 
shifted arbitrarily along the vertical axis are shown. The upper 
curve is an integral light curve of the four quasar images (see 
the text for more details). 



to note as well, that the errors indicated in Table ^reveal no 
seeing dependence for PSFs up to 1.2". However, one can eas- 
ily notice much larger errors for three dates, when the seeing 
was 1.4" to 1.5". 

The Q2237+0305 A,B,C,D light curves are shown in FigQ] 
To better illustrate their similarity, they are arbitrarily shifted 
along the magnitude axis. The upper curve is a total light curve 
obtained from the aperture photometry of the central part of 
the object within a diaphragm of 3"radius, and with the subse- 
quent subtraction of the intensity contribution from the galaxy. 
Polynomial fits (fifth order polynomials were used) are also 
shown. All light curves are well correlated, thus demonstrating 
extremely low microlensing activity during this time period, 
while the quasar revealed rather high variability, - 0.15 m be- 
tween the Julian dates 860 and 910, i.e. about 0.003'" per day. 
This is almost five times more rapid brightness change as com- 
pared to the measurements by 0stensen et al. ( 1996 1 in 1994. 
Similar synchronous brightness variations can be seen in the 
V light curves obtained by the OGLE group for the same time 
period, (http://www.astrouw.edu.pl ogle/ogle3/huchra.htmli, 
though the authors do not note similarity between the light 
curves of the components. A quantitative measure of similarity 
between the light curves can be seen from Table where 
the corresponding correlation indices are shown for our 
photometry and for the OGLE data. Since the expected time 
delays are very short, the table justifies an attempt to undertake 
calculation of the time delays, which is possible thanks to 
subdued, though non-zero microlensing activity. 

It is interesting to estimate the level of possible microlens- 
ing activity during this time period. In Fig|2] variations of mag- 
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Table 1. Photometry of Q2237+0305 A,B,C,D in the R band from observations with the Maidanak 1.5-m telescope in June- 
October 2003. 



Date 



JD 

(2452000+) 



B 



D 



FWHM 

(arcsec) 



Number 

of 
frames 



jun27 
jun28 
jun29 
julOl 
jul04 
jul07 
jul08 
jullO 
jull4 
jull9 
jul26 
jul28 
jul30 
jul31 
aug05 
aug06 
aug08 
aug09 
auglO 
augll 
augl4 
augl5 
augl6 
aug21 
aug22 
aug23 
aug24 
aug25 
aug26 
aug27 
aug28 
sepOl 
sep02 
sep05 
sep06 
sep 12 
sep 13 
sep 15 
sep 16 
sep 17 
sep 18 
sep 19 
sep20 
sep21 
sep24 
sep25 
sep27 
sep29 
sep30 
octOl 
oct02 
octll 
octl2 



~8T8~ 
819 

820 
822 
825 
828 
829 
831 
835 
840 
847 
849 
851 
852 
857 
858 
860 
861 
862 
863 
866 
867 
868 
873 
874 
875 
876 
877 
878 
879 
880 
884 
885 
888 
889 
895 
896 
898 
899 
900 
901 
902 
903 
904 
907 
908 
910 
912 
913 
914 
915 
924 
925 



17.055+0.003 
17.071+0.014 
17.069+0.013 
17.076+0.005 
17.049+0.008 
17.053+0.010 
17.070+0.009 
17.052+0.014 
17.077+0.005 
17.086+0.006 
17.077+0.006 
17.057+0.018 
17.056+0.003 
17.069+0.001 
17.073+0.007 
17.070+0.006 
17.070+0.004 
17.072+0.005 
17.069+0.005 
17.066+0.018 
16.925+0.009 
17.019+0.006 
17.089+0.006 
17.062+0.010 
17.032+0.004 
17.048+0.005 
17.038+0.009 
17.037+0.005 
17.017+0.006 
17.039+0.007 
17.026+0.005 
17.006+0.010 
16.978+0.032 
16.981+0.012 
16.996+0.016 
16.826+0.015 
16.991+0.003 
16.924+0.024 
16.923+0.006 
16.975+0.003 
16.953+0.008 
16.946+0.004 
16.937+0.003 
16.923+0.005 
16.924+0.006 
16.933+0.004 
16.928+0.008 
16.917+0.007 
16.927+0.002 
16.932+0.005 
16.929+0.003 
16.938+0.003 
16.920+0.005 



18.549+0.029 
18.615+0.030 
18.593+0.026 
18.604+0.023 
18.533+0.032 
18.574+0.004 
18.590+0.054 
18.598+0.035 
18.605+0.025 
18.595+0.029 
18.594+0.011 
18.622+0.015 
18.585+0.005 
18.566+0.012 
18.564+0.016 
18.592+0.009 
18.570+0.007 
18.578+0.013 
18.577+0.017 
18.606+0.046 
18.454+0.059 
18.534+0.026 
18.622+0.018 
18.543+0.015 
18.555+0.012 
18.567+0.005 
18.533+0.015 
18.546+0.012 
18.525+0.006 
18.555+0.023 
18.541+0.019 
18.535+0.022 
18.583+0.026 
18.515+0.029 
18.491+0.065 
18.567+0.227 
18.533+0.031 
18.450+0.039 
18.513+0.028 
18.471+0.007 
18.499+0.010 
18.473+0.010 
18.466+0.004 
18.465+0.010 
18.487+0.011 
18.470+0.013 
18.503+0.019 
18.440+0.010 
18.469+0.003 
18.470+0.011 
18.476+0.008 
18.479+0.011 
18.433+0.009 



18.285+0.013 
18.282+0.012 
18.286+0.015 
18.272+0.021 
18.287+0.025 
18.260+0.005 
18.274+0.032 
18.310+0.034 
18.305+0.023 
18.280+0.023 
18.289+0.020 
18.280+0.057 
18.269+0.022 
18.284+0.015 
18.262+0.012 
18.246+0.014 
18.259+0.008 
18.270+0.016 
18.276+0.015 
18.236+0.036 
18.201+0.022 
18.235+0.030 
18.288+0.012 
18.242+0.017 
18.245+0.007 
18.249+0.012 
18.246+0.006 
18.237+0.012 
18.222+0.007 
18.258+0.010 
18.222+0.004 
18.193+0.017 
18.183+0.059 
18.219+0.024 
18.169+0.057 
18.089+0.023 
18.215+0.008 
18.103+0.046 
18.154+0.033 
18.201+0.004 
18.157+0.011 
18.142+0.007 
18.151+0.012 
18.123+0.016 
18.144+0.003 
18.135+0.010 
18.139+0.025 
18.115+0.010 
18.133+0.007 
18.138+0.015 
18.126+0.006 
18.137+0.009 
18.117+0.007 



18.392+0.031 
18.389+0.014 
18.450+0.021 
18.437+0.042 
18.413+0.042 
18.415+0.008 
18.424+0.109 
18.414+0.036 
18.422+0.045 
18.436+0.021 
18.439+0.008 
18.446+0.049 
18.430+0.034 
18.413+0.012 
18.416+0.016 
18.425+0.012 
18.413+0.012 
18.437+0.014 
18.435+0.020 
18.440+0.043 
18.437+0.039 
18.362+0.019 
18.427+0.022 
18.385+0.018 
18.399+0.015 
18.400+0.019 
18.413+0.014 
18.401+0.022 
18.380+0.017 
18.388+0.011 
18.394+0.011 
18.390+0.009 
18.373+0.047 
18.366+0.036 
18.282+0.019 
18.394+0.153 
18.327+0.033 
18.340+0.051 
18.395+0.048 
18.310+0.014 
18.358+0.014 
18.319+0.012 
18.305+0.013 
18.308+0.011 
18.325+0.011 
18.331+0.010 
18.311+0.016 
18.329+0.010 
18.313+0.004 
18.323+0.006 
18.329+0.004 
18.315+0.009 
18.287+0.012 



~nr 

1.0 
1.0 
0.9 
1.0 
0.9 
1.0 
1.2 
1.1 
0.9 
1.0 
1.2 
1.1 
0.9 
1.2 
0.9 
0.9 
0.8 
0.9 
1.0 
0.9 
1.1 
1.1 
0.8 
0.9 
1.0 
1.0 
1.0 
0.8 
0.8 
1.0 
1.0 
1.5 
1.1 
1.1 
1.4 
0.9 
1.5 
1.2 
0.9 
0.9 
1.1 
0.8 
1.2 
1.0 
0.9 
1.2 
1.2 
0.8 
0.9 
0.8 
0.9 
1.0 



6 
4 
4 
4 
6 

12 
4 
9 
6 
4 
5 
4 
4 
4 
4 
4 
6 
4 
4 
4 
4 
6 
4 
4 
4 
4 
4 
3 
8 
4 
12 



nitudes of individual quasar components with respect to their 
integral light curve are shown, which are probably caused by 
microlensing. The figure illustrates, that the microlensing vari- 
ations are rather slow and may be regarded as linear trends 
during the 100-day monitoring period. The mean slopes fit- 
ted to these microlensing residuals were -0.000136, 0.000125, 
-0.000251, and 0.000264 magnitudes per day for the A, B, C 
and D components, respectively. These are important quanti- 
ties, which we will need in the next section to analyse uncer- 
tainties of the time delays estimates. 



We tried to fit the observed brightness changes by polyno- 
mials of the 3-rd, 4-th, 5-th, and 6-th orders, and found that 
there is no justification to use the polynomial degree higher 
than four in calculations, since scatter of observational data 
points from the fitting curves decreases very slowly as the 
polynomial order grows. For example, the fourth-order polyno- 
mial gives the RMS deviations for the datapoints of 0.0148 m , 
0.0228 m , 0.0161 m , 0.0206 m , and 0.0098 m for the A, B, C 
and D components and the quasar, respectively, whereas with 
the fifth-order polynomial, the RMS deviations are 0.0147" 1 , 
0.0227'", 0.0159 m , 0.0204'", and 0.0098 m . 
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Table 2. Correlation indices between the light curves of the 
Q2237 components calculated for our data, and for the OGLE 
photometry. Quantities in the last line and column of the first 
table, both marked as QSO, are correlation indices between a 
particular light curve and the integral light curve of the quasar. 

Our data 



Component 


A 


B 


C 


D 


QSO 


A 


1 


0.899 


0.960 


0.874 


0.978 


B 


0.898 


1 


0.888 


0.887 


0.924 


C 


0.960 


0.888 


1 


0.874 


0.968 


D 


0.874 


0.887 


0.874 


1 


0.922 


QSO 


0.978 


0.924 


0.968 


0.922 


1 


OGLE data 


Component 


A 


B 


C 


D 




A 


1 


0.945 


0.953 


0.798 




B 


0.945 


1 


0.914 


0.754 




C 


0.953 


0.914 


1 


0.726 




D 


0.798 


0.754 


0.726 


1 
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Fig. 2. Residual variations of brightnesses of the A-D images 
taken from the light curves of June 27 - October 12, 2003 
(Fig.l) by subtracting the integral quasar light curve. These 
residuals may be treated as a result of microlensing. 



This quantities may be regarded as an upper estimate for 
the actual errors of our photometry, while the values 0.008" 1 , 
0.0 18 m , 0.016'", and 0.020'" indicated in the second paragraph 
of this section, should be regarded as a measure of the intrin- 
sic convergence of the photometry for images A, B, C and D, 
respectively. Thus the data presented in Table ^ an d shown in 
Figffl are accurate enough to attempt calculations of time de- 
lays. Table|2]contains an additional argument in favour of such 
an attempt, - the light curves of all the four quasar images are 
well correlated. 

Nevertheless, to be sure that such featureless light curves 
are really capable of providing reliable estimates of time de- 



lays, to investigate the effects of photometry errors and residual 
microlensing effects, and to obtain consistent estimates of the 
time delay uncertainties, numerical simulations involving over 
2000 statistical trials were undertaken. The simulation results 
will be briefly described at the end of the next section, after the 
methods used to determine time delays are presented. 

3. Time Delay Estimation 

The detection of correlated variations of brightnesses of the 
Q2237 components is important in itself, since it is the most 
convincing evidence of the gravitational lensing nature of the 
system. The deeper value of these observations is in the possi- 
bility to estimate the time delays between the brightness fluc- 
tuations seen in different lensed quasar images. 

Estimation of time delays from the actual light curves, with 
their inevitable uneven sampling and poor statistics, is a rather 
difficult problem. A great variety of approaches has been elab- 
orated, with the most important ones being analysed by Gil- 
Merino et al. (2002 1. All of these algorithms are devised for 
the situation when the characteristic timescale of the quasar 
brightness variations and sampling time intervals of observa- 
tions are much less than the expected time delays. The existing 
macrolensing models for this system predict time delays which 
do not exceed a day. Therefore, we have a somewhat atypi- 
cal situation, - we are trying to estimate the time delays from 
brightness variations with the characteristic timescale much 
larger than the expected value of the time delay. Moreover, 
we used photometry data averaged within a night, and thus the 
sampling time interval exceeds the expected values of time de- 
lays too. 

Taking all this into account, we used the following ap- 
proach. We admit that the quasar brightness variations are 
rather slow and can be represented accurately enough by a 
fourth or fifth order polynomial. Then, representing bright- 
ness variations of a component, e.g. A, with a smooth curve 
we may calculate a normalized cross-correlation func- 
tion (NCCF) of this smooth curve with the actual light curve of 
another component, say, B, represented by discrete data points: 



k{r) = 



jgn \Mi + t) - f A (T)][m B (ti) - m 

(vy m y> 2 



(i) 



where are discrete brightness samples at the time mo- 

ments tj, f&(ti + t) are corresponding values of the approximat- 
ing polynomial at the time moments f; + r, f A and JEb are cor- 
responding average values, and t is the time shift. The average 
values /(r), m and variances <r^,cr^ are calculated according to 
formulas: 



i - i " 

- X fiU + r), m = - V m(ti); 



•f 



w-Z 



[/(f,+r)-/(r)] 2 
n-l 



[m(tj) — m] 2 



(2) 



(3) 



With such an approach, we expect to find cross-correlation 
peaks at the cosmological lags. We adopted this procedure be- 
cause it is less sensitive to bias by any false peaks for lag 
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that can result from the kinds of correlated brightness fluctua- 
tions at lag when working with discrete data available from 
optical monitoring. These correlated photometry errors are ex- 
pected from the inevitable residual sensitivity of the photome- 
try to seeing effects, that is often observed, (e.g. Vakulik et al. 
2004 1. The effect of such "zero lag" is also discussed in Colley 
et al.(2003) and interpreted as a "frame-to-frame correlation 
error in the photometry". 

The NCCF, calculated in this way, is a smooth function, 
with one maximum, which changes very slowly, and does not 
have other local extrema at least within -10 < r < 10 days. It 
allows us to unambiguously interpret the location of the maxi- 
mum of this function as an optimal value of the lag, which we 
accept as the time delay between the brightness variations of 
two components. In this way we can estimate the time delays 
for the light curves of all the components with respect to the 
fitting curve of the A component. 

Similar calculations can be made with respect to the fitting 
curves of other components, and thus a complete set of the time 
delay estimates will be obtained. Such estimates cannot be re- 
garded as independent ones however, since they are affected 
with the same random errors of light curves of the "template" 
components, while their differences are caused by the differ- 
ence between the approximating curves used. 

We also propose a more general approach, consisting in a 
joint fitting of all the four light curves by a single polynomial, 
with a normalizing factor, magnitude shift, and time shift of 
each particular light curve being treated as unknown parame- 
ters and found from a condition: 

N 

f(a,ti,Tj) = Y j an{t i + T j f. (5) 

n=l 

Here, i is a number of a data point in the light curve of 
the j-th component, kj is a normalizing factor for the j-th light 
curve, and tj is the unknown time delay for the j-th component. 
The function f(a, f,, Tj) is an A^-th order polynomial without 
the zeroth order term which is represented for each component 
separately as the magnitude shift aoj. 

These two approaches, - which we marked as NCCF and 
JFM (joint fitting method), - have been applied to the photome- 
try data of June-October 2003, presented in Table 5 and shown 
in Fig.^ The estimates of time delays of the B, C and D images 
with respect to image A are presented in Table [3] Column 2 
contains the estimates, taken with the NCCF method. They are 
the results of averaging the time delay estimates calculated with 
respect to the fitting curves of all the components separately. 
In column 3, results of applying the JF method are shown. A 
comparison of these two columns illustrates the validity of both 
methods, - they produce almost the same time delay values. 

If, in addition to the quasar intrinsic brightness changes, 
microlensing brightness variations took place during our moni- 
toring time period, this might cause certain bias in the time de- 
lays estimates, and difference in the amplitudes kj determined 
from (4). With the assumption that there were no strong mi- 
crolensing events during our monitoring, possible minor effects 



of microlensing at the time interval of 100 days can be approx- 
imately represented by linear trends (see Fig0. To check a 
suggestion that the observed amplitude difference results from 
microlensing of the components, we changed Eq. (4), having 
introduced, in addition to the constant terms aoj, linear terms 
of the polynomial as well, a\j, for each component separately: 

2 th T >) + a °-j + fll J^» + T >) ~ m / f ')] 2 = min - ( 6 ) 

With this modification of the JFM, which we called JFMi, 
the averaged linear terms of the components would represent a 
linear term of the quasar brightness variation, while their differ- 
ence might be interpreted as the effects of weak microlensings 
of individual components. If this were a case, such a procedure 
would equalize the amplitudes of the components, but this does 
not happen. The amplitudes can also be equalized forcedly by 
assuming kj = 1 for all the j-s in condition (6), (JFM2 modifi- 
cation of the JFM). After such a procedure, however, residuals 
in Eq.4 became much larger than for the JFM and JFMi, - the 
light curves just stopped to be similar, - and this finally rejects 
a possibility for the amplitudes to be equalized by introducing 
microlensing trends. 

Therefore, our assumption that the quasar brightness vari- 
ation is observed in different components with different am- 
plitudes because of microlensings does not get a solid confir- 
mation. It should be note also, that the time delays calculated 
with the JFMi are less reliable as compared to those obtained 
with the JFM, since the JFMi deals with a larger number of un- 
known parameters. Taking all this into account, we placed the 
results of applying the JFMi and JFM2 in Table [3] (columns 4 
and 5), but did not use them in calculations of the average time 
delays (the bottom of Table [3}. Possible reasons for the ob- 
served difference in light curves amplitudes will be discussed 
in Section 5. 

A 4-th order polynomials were used for fitting in all cases. 
Further increasing of the polynomial order does not improve 
the accuracy of fitting, and further simulation shows that the 
random error of time delay determination increases slightly as 
the polynomial degree grows. 

To determine reliability of estimates, and to analyse the 
effects of random photometric errors, numerical simulations 
were undertaken. The eighth-order polynomial fit to the im- 
age A light curve was adopted as an intrinsic quasar brightness 
curve template. Thus we admitted, that the actual brightness 
variation of the quasar might be more complicated as com- 
pared to the 4-th order polynomial fit used for time delay es- 
timations. Individual light curves of the components were then 
constructed by adding to the template curve, at the time mo- 
ments corresponding to those of observations, the random noise 
quantities with the RMS values of 0.015 m , 0.023" 1 , 0.016 m , and 
0.021'" inherent in photometry for images A, B, C and D, re- 
spectively. A zero time lag was adopted for all the images, 
and 2000 random samples for each quaternary of images were 
formed for the simulation. The algorithms described in the pre- 
vious section were applied to the random samples to calculate 
the average time delays with respect to image A, and their RMS 
deviations. The estimates of time delays from these simulated 
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Table 3. Estimates of time delays for Q2237+0305 system 
with respect to the A image, (hours), taken with different ap- 
proaches from the Maidanak data, and from those of the OGLE 
group. The estimates averaged over NCCF and JFM for both 
Maidanak and OGLE data are also shown, (weighted average 
values). The JFMi and JFM2 columns are shown just to illus- 
trate the attempt to equalize the amplitudes of the light curves, 
and they were not used to calculate the average time delays. 



Maidanak data 



Im. pairs 


NCCF 


JFM 


JMF, 


JFM, 


BA 


-13.2 


±59 


-16.1 ±62 


33.1 


-36.2 


CA 


-31.9 


±42 


-37.2 ±42 


36.9 


-10.6 


DA 


6.0 i 


:60 


5.0 ± 63 


9.1 


-40.1 


OGLE data 


BA 


1.6 i 


:52 


5.0 ± 57 




7.2 


CA 


105.8 


±40 


105.1 ±43 




117.8 


DA 


-5.9 d 


= 152 


-3.6 ± 156 




17.8 


Average 












BA 


-6 + 


:41 








CA 


35 ± 


30 








DA 


2± 


44 









Table 4. Model predictions for the time delays tba, tca, and 
tda for Q2237+0305 system, (hours). 



Reference 


Lens model 




TCA 


Tda 


Schneider, 


Constant mass-to-light 








etal.(1988) 




-2.4 


-29.5 


-26.6 


ratio 








Rix et al. 


Model 1 


2.1 


11 


7.1 










(1992) 


Model 2 


1.5 


10.1 


6.1 












Model 2a 


-1.7 


9.8 


3.7 


Wambsganss 


Point lens 


-2.97 


17.41 


4.87 


& Paczynski 


Isothermal sphere 


-1.51 


8.91 


2.46 


(1994) 


Best fit 


-0.44 


2.54 


0.7 


Schmidt, 


Bar accounted 


-2.0 


16.2 


4.9 


etal.(1996) 











samples did not reveal any noticeable biases, and the RMS er- 
rors turned out to be about 62, 42 and 63 hours for B, C and D, 
respectively. 

We have also considered the possibility that a small amount 
of systematic microlensing could influence these time delay 
determinations. We have simulated the effects of such small 
systematic microlensing as a sustained brightness increase (de- 
crease) during 100 days of our monitoring. To do this, the ac- 
tual light curves were distorted by linear trends, and the result- 
ing time delay biases were estimated. For our low amplitude 
trends, the biases of the time delay estimates were found to de- 
pend almost linearly on the value of trend, with 1 day for the 
daily rate of 0.0005"'. In Sec. 2, (see also Fig. 2) the maximal 
possible value of a microlensing trend has been estimated to 



be about 0.00025™ per day. Therefore, the maximum possible 
time delay bias caused by uncompensated microlensing events 
does not exceed 12 hours. 

We also tried to calculate the time delays using the OGLE 
data, (http://www.astrouw.edu.pl ogle/ogle3/huchra.htmli. The 
results of applying the NCCF and JF methods are shown in 
Table |3 (columns 2 and 3), as well as the attempts to equal- 
ize the amplitudes of the OGLE light curves, (JFM2 column). 
The results show, that the high accuracy inherent the OGLE 
photometry, does not seem to improve the results. The error es- 
timates resulting from our simulations as above, are still very 
large, and the time delays differ from our values and also from 
the values predicted from available macro-imaging models. 

The results of averaging the estimates taken with the NCCF 
and JFM for both the Maidanak and OGLE light curves are also 
shown in Table|5] where the weighted average values are given. 
In the following sections, we analyse these results and discuss 
the pitfalls and prospects of observational measurement of time 
delays for the Q2237+0305 system, from daily sampled optical 
wavelength data. 

4. Time delay predictions for Q2237 from 
macro-imaging models 

The first predictions of time delays for the Q2237+0305 sys- 
tem appeared soon after its identification as a quasar quadru- 
ply imaged by the gravitational field of a foreground galaxy, 
(Schneider et al. 113551 Kent & Falco [13551 . Schneider et al. 
(1988 1 treated the lensing galaxy projected mass distribution 
as proportional to the observed light distribution and found that 
"a constant mass-to-light ratio, elliptical, de Vaucoulers bulge 
reproduces the observed configuration remarkably well." The 
model was improved in 1992 with the use of a more accu- 
rate galaxy light distribution using the Hubble Space telescope 
data (Rix et al. I1992> . Kent and Falco ([1988) used another 
approach and considered elliptical mass distributions with a 
de Vaucouleurs and a King profile mass distribution. Another 
way to account for the mass distribution ellipticity was pro- 
posed, e.g., by Kochanek J1991I and Wambsganss & Paczynski 
(1994 1. They used a shear, y, while the mass distribution in 
the lens was adopted to be spherically symmetric, - a singular 
isothermal sphere and a point mass plus different types of the 
shear, (Kochanek 1991 1, and a singular power-law sphere with 
an external shear (Wambsganss & Paczynski 1994 1. A massive 
disk with a central nucleus as a point mass, isothermal sphere 
and Plammer's sphere were also considered as models for the 
Q2237 lensing galaxy, (Minakov & Shalvanin ll991> . An at- 
tempt to account for the galaxy bar should be mentioned also, 
(Schmidt et al. 119981 . They used two-dimensional Ferrers pro- 
files for the mass distribution in the bar, and a power-law ellip- 
tical mass distribution for the bulge, with the bar representing a 
small perturbation of the deflecting field of the bulge and con- 
stituting about 5% of the bulge mass in the region inside the 
quasar images. 

Predictions of the time delays for the Q2237 images that we 
could find in the works cited above, are collected into Table |4] 
with the proper referencing and brief comments concerning 
each particular model used. Here, Models 1 and 2 by Rix et 
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Fig. 3. light curve of image A from observations on 
Maidanak in June 27 - October 12, 2003, (our photometry), - 
upper curve, solid circles, - and V light curve of image A built 
from the OGLE data, - bottom curve. 

al.(|T992i mean, respectively, R 1 ^ 4 profile, R 1 ^ 4 plus unresolved 
nucleus, and Model 2a means Model 2 with only image posi- 
tions matched. The predicted values of time delays are seen to 
be obviously model dependent. The values of the expected time 
delays calculated by Wambsganss & Paczynski ( 1994 1 exhibit 
a tendency to decrease as the index of power f3 in their power- 
law mass profile grows - from f3 — for a "point lens", through 
y8 = 1 for an isothermal sphere, and up to j3 = 1.71, ("best 
fit model"). The B component is always leading, in contrast to 
the model by Rix et al. (I1992> . A rather good agreement of 
their point lens model predictions with those by Schmidt et al. 
( 11998 1 should be noted, where the effect of the galaxy bar was 
taken into account. 

In light of the errors and simulations presented in section 3, 
we can make no conclusions about the comparison of models 
with our lag calculations from monitoring data. 

5. Quasar Intrinsic Brightness Fluctuation 
Amplitude 

In analysing our calculations, one more thing should be noted 
and discussed. It has been shown in Sec. 2 (see Table[2}, that the 
light curves of all the components are highly correlated, both in 
the R filter (our observations), and in the V band (OGLE data). 
It should be noted however that the R light curves built from 
the same Maidanak data by the Moscow group, (Koptelova et 
al. 2005 1 differ from our light curves, with the cross correlation 
indices being much lower than those shown in Table [2] and 
ranging from 0.2 to 0.7. 

If the light curves observed during the time period from 
June 27 to October 12, 2003, represent a quasar intrinsic bright- 
ness change, then its amplitude seen in different lensed images 
is believed to be the same. As seen from Table [5] where the 
amplitudes of our R light curves and those of the OGLE group 
in the V band are presented in the units of "our" A component 



Table 5. Relative amplitudes of R light curves, (our data), and 
light curves in V (OGLE data) calculated for the time period 
June 27 - October 12, 2003. 



Image 


R 


V 


A 


1.000 


1.263 


B 


0.809 


0.796 


C 


1.004 


0.854 


D 


0.751 


0.370 



amplitude, this is not the case. The amplitudes of the A and C 
light curves are noticeably larger than that of D both from our 
data and from the OGLE group data. Inconsistencies between 
our measurements and the OGLE data can also be seen from 
this table. To better demonstrate what kind of inconsistencies 
between the Maidanak and OGLE data take place, - and, at the 
same time, how similarly they look, - we plotted our R and the 
OGLE V light curves of image A in the same picture, (Fig[3}. 

A possible explanation of these complex differences in fluc- 
tuation amplitude may be the following. We adopt the pic- 
ture that all four image components may be microlensed but 
far enough from a caustic in the microlensing pattern that al- 
most no change in relative brightnesses of the 4 images is seen 
during our monitoring period. It is also possible that the mi- 
crolensing caustic pattern in the four images is moving ap- 
proximately parallel to the local caustic structure. It is known 
that microlensing can cause an image to become bluer as the 
caustics cross the quasar structure and magnify most strongly 
the bluest innermost region (Wambsganss & Paczynski. [T99T1 
Vakulik et al., 2004 1. Presumably at the 2003 time of these ob- 
servations, image A was brightest (and image B was faintest) 
because of microlensing image magnification (demagnifica- 
tion). Such magnification can alter the color of the images as 
noted previously, and depending upon the details of which part 
of the image are most strongly microlensed, and which part of 
the quasar became brighter in the course of the quasar intrin- 
sic brightness fluctuation, color microlensing effects of several 
kinds might be expected. These topics have not been explored 
with extensive simulations to date, because there is as yet no 
standard quasar model. We also note the possibility that the un- 
equal amplitudes might result in part from complex structure 
in either the quasar image (Schild & Vakulik 2003 1, or in the 
lens galaxy. In the latter case, a subtle unresolved lens galaxy 
brightness feature might cause the erroneous estimation of an 
individual image brightness; however we would expect this to 
have been noticed in other Q2237 data analysis. 

We also note a possible artifact from systematic errors in 
image processing have been discussed already in a number of 
publications, e.g., Alcalde et al. (2002i, Burud et al. ( 1998 »; 
Vakulik et al. (2004i. Because of closeness of quasar images to 
the galaxy nucleus and to each other, errors of processing due 
to incorrect subtraction of foreground galaxy light or due to in- 
correct subtraction of light from the other images are possible, 
different for different photometry methods. 
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6. Strategy to Detect Time Delays 

The interesting X-ray observation of possible 2.7 hour time de- 
lay between the A and B images (Dai et al. 2003) raises the 
question of limits of our ability to measure with precision the 
lags for all four images from optical data. In spite of years 
of Q2237 brightness monitoring, little is yet known about the 
quasar's intrinsic brightness fluctuations on long or short time 
scales. 

On short time scales, Burud et al. (1998 1 and Cumming 
et al. dl995l monitored Q2237 intensively for three nights 
with the NOT and with the CFHT. And nightly monitoring by 
Vakulik et al. (120041 for 46 nights in 2000 with the R filter 
precludes fluctuations of more than 1%. 

We can also recall what is known about rapid fluctuations 
in other quasars. Colley & Schild (2003) present convincing 
evidence for Q0957 nightly brightness fluctuations at the level 
of 0.01 magnitude (their Fig. 1). Gopal-Rrishna et al. (|2003l 
report detection of intranight variations of brightness of about 

0. 01 - 0.02 magnitudes in their sample of radio-quiet quasars. 
At the same time, in a survey of 23 high-luminosity (-27 < 
M v < -30) radio-quiet quasars by Rabbette et al. {1998 1, "no 
evidence for short-term variability greater than about 0. 1 mag- 
nitudes was detected in any of the 23 sources". 

Long-term variability of quasars is rather well documented 
with the use of huge amounts of observational data, - e.g., re- 
sults of statistical processing of a sample containing over 40000 
quasars as reported by de Vries et al. ( 2005 1. Structure function 
analysis of this sample produced estimates of short-term vari- 
ability consistent with the above, and also produced some im- 
portant statistical characteristics of the variability of quasars, 
with the following ones relevant here: l)the increase of vari- 
ability with increasing time lags is monotonic and constant up 
to time-scales of ~ 40 years; 2)variability increases toward the 
blue part of the spectrum, and 3)high-luminosity quasars vary 
less than low-luminosity quasars. 

Using this information, and based upon the results of this 
paper, we propose the following strategy for future attempts to 
observationally determine the time delays for the Q2237+0305 
gravitationally lensed quasar: 

- since variability of quasars tends to increase towards the 
shorter wavelength, that is found for a sample of quasars by 
de Vries et al. (2005 1, and also detected for the Q2237 sys- 
tem, (Vakulik et al. 2004 1, observations should be shifted 
towards as much shorter wavelengths as possible; 

- monitoring in sets with durations of several weeks and sam- 
pling interval of hours; 

- observations should be made under the image quality of 
0.6-0.7"and better, and with the photometry accuracy of 
0.005'" and higher. 

7. Conclusions 

1 . Variations of brightness of the four quasar images of Q2237, 
observed during almost 100 days in June - October 2003, 
are due mainly to the quasar intrinsic brightness change. 
Microlensing activity was substantially subdued during this 
time period and did not exceed 0.00025 magnitudes per day. 



2. The estimates of time delays obtained with the two methods 
proposed, are consistent with each other, but differ from those 
predicted by the existing macro-imaging models. The resulting 
error bars neither confirm nor rule out any of the competing 
models. 

3. The ultimate limits of time delays estimation from the avail- 
able observations were tested in simulations and have shown 
that, with the particular type of the quasar variability, particular 
photometric errors and contribution from microlensing activity, 
one should not expect better results and higher accuracy. 

4. The measured time delays between the 4 quasar images 
are much shorter than the time scales of observed optical mi- 
crolensing brightness changes. Therefore the long history of 
analysing brightness changes in the Q2237 literature is based 
upon a correct belief that the observed fluctuations are due 
mainly to microlensing events. However, the detected bright- 
ness change of the quasar, - 0.15'" between the Julian dates 
860 and 910, i.e. about 0.003'" per day, - demonstrates that the 
quasar intrinsic variability may noticeably contribute to the ob- 
served light curves of the components, and thus it should be 
taken into account when analysing and interpreting statistics of 
microlensing in the system. 
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